THE ATIYAH-HITCHIN BRACKET FOR THE CUBIC 
NONLINEAR SCHRODINGER EQUATION. 

I. GENERAL POTENTIALS. 



K.L. VANINSKY 

Abstract. This is the first in a series of papers on Poisson formalism for the 
cubic nonhnear Schrodinger equation with repulsive nonlinearity and its relation 
to complex geometry. In this paper we study general continuous potentials. 
We demonstrate that the Weyl functions of the corresponding auxiliary Dirac 
spectral problem carry a natural Poisson structure. We call it the Atiyah- 
Hitchin Poisson bracket. We show that the Poisson bracket on the phase space 
is the image of the Atiyah-Hitchin bracket on Weyl functions under the inverse 
spectral transform. 



1. Introduction. 

1.1. Statement of the problem. All 1+1 differential equations like the Korteweg- 
de-Vries, the modified Korteweg-de-Vries, the cubic nonlinear Schrodinger equa- 
tion, the Toda lattice, the Camassa-Holm equation etc., which are analyzed by 
the inverse spectral transform, are Hamiltonian systems. 
The cubic NLS with repulsive nonlinearity^ 

i^' = _^" + 2|V'|V, 

where ip{x, t) is a complex function, will serve as our model example. We consider 
this problem on the entire line, i.e., x G M}. We assume that the phase space 
A4 consists of functions ip{x) and we do not impose any condition on ip{x) except 
continuity. We restrict our attention to some specific function classes, e.g., periodic 
or rapidly decaying, in subsequent papers. To make apparent the algebraic nature 
of our considerations, we assume that all functionals A, B : M. —>■ C are Frechet 
differentiable, all integrals converge, etc. 

The cubic NLS equation is a Hamiltonian system 

r = {^,n}, 

with the classical bracket 

{A,B} = 2t^ — — (1-1) 

J 5tp{x) oil){x) dil}{x) 6tp{x) 



The work is partially supported by NSF grant DMS-9971834. 

^Prime ' signifies the derivative in the variable x and dot • the derivative with respect to time. 
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and Hamiltonian 



H 



+ 



^ dx. 



The NLS equation arises as a compatibility condition for the commutator re- 
lation for some specially chosen differential operators. This leads to an auxiliary 
linear spectral problem for the Dirac operator 



1 
-1 





itp 



acting in a space of vector functions 



fi{x,X) 

f2{x,X) 



The goal of this and subsequent papers is to demonstrate that the Hamiltonian 
formalism (Poisson bracket) is built into the complex geometry of the Dirac spectral 
problem. 



1.2. Description of results. We associate to a potential il){x) of the Dirac spec- 
tral problem a pair (F, X), i.e., we have a map: 

M {t,a:). 

We call this map the direct spectral transform. In this pair F is a two sheeted 
covering of the complex plane of the spectral parameter cut along the real line 
(see Figure 1). For any point Q e F, its projection to the plane of the spectral 




parameter A is denoted by X{Q). The holomorphic function X{Q) defined on F is 
the so-called Weyl function of the Dirac spectral problem. We show that the pair 
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(r, X) carries a natural Poisson structure. We call it the Atiyah-Hitchin bracket. 
Apart from insignificant constants the AH bracket is given by the formula 

This is the first theorem of the paper. 

It turns out that the direct spectral transform can be inverted 

iT,X) M. 

We call this map the inverse spectral transform. We show that the image of the 
AH bracket under the inverse spectral transform becomes the bracket 1.1 on phase 
space. This is our second theorem. 

The pair (F, X) introduced in this paper is a skeleton of an analogous construc- 
tion for periodic and scattering potentials in subsequent papers [14] and [15]. In 
these cases there is an additional structure which allows one to identify the values 
of the meromorphic function X on different banks of the cut. 

The question of construction of canonical coordinates can not be resolved with- 
out assumptions on a potential. Canonical coordinates are given in terms of sin- 
gularities of the Weyl functions on the real line, or in other words in terms of the 
spectrum of the auxiliary operator. At the moment the only equations which have 
been analyzed are the finite Toda lattice, [11], and the Camassa-Holm equation 
with rapidly decaying initial data of one sign, [12]. For both of these systems the 
spectrum of the corresponding auxiliary spectral problem is an isolated discrete 
set^ . We refer to [11] and [12] for details. 

1.3. Historical remark. Another approach which relates Hamiltonian theory and 

complex geometry was developed recently. 

In a remarkable paper Krichever and Phong [6] proposed a construction of a 
symplectic formalism for integrable equations with periodic initial data. For the 
latest exposition of their results and connections with Seibcrg-Witten theory see [2]. 
The Krichever- Phong approach is designed for 2+1 systems like the Kadomtsev- 
Petviashvili equation and the 2D- Toda lattice. It can also be extended to the 
scattering case, [13]. In the simplest cases the simplectic form can be inverted ex- 
phcitly to obtain the Poisson bracket. For example the Gardner-Zakharov-Faddeev 
bracket for KdV or classical bracket 1.1 (both are constant on the phase space) can 
be obtained this way. In the case of a bracket with variable coefficients the situation 
is different. It is not possible to invert explicitly the corresponding higher (second) 
symplectic form even for the case of NLS with periodic boundary conditions, see 
[13]. These prompt us to find another approach to the Poisson formalism. 



For the Camassa-Holm equation and the open Toda lattice the corresponding Poisson bracket 
on the Weyl function is also given by formula 1.2 or its' reduction. 
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1.4. Content of the paper. In Section 2 we define tiic Atiyah-Hitchin bracket 
for rational maps of CP^ into itself. We state its basic properties and prove that 
the AH bracket is invariant under linear fractional transformations. The direct 

spectral transform and inverse spectral transform are defined in Section 3. We 
study the Poisson bracket in Section 4, which consists of two parts. In the first 
part we compute the image of the classical bracket for the Weyl function. In the 
second part we show that the bracket on phase space is an image of the AH bracket 
under the inverse spectral transform. 

2. The Atiyah-Hitchin bracket. 

Let X{X) : CP^ CP^ be the rational function of degree N such that 

• A'(oo) = 0, 

• the preimage X~^{oo) consists of exactly N points. 

To any A'"-monopole solution mjv(x), x e M^, of the Bogomolny equation Atiyah 
and Hitchin, see [1], associate the scattering function X(X), with these properties. 
Thus we have the direct spectral transform 

This map is injective due to a theorem of S. Donaldson, [1]. 
The function A* (A) has the form 



N-l 
1=1 



N-l 

A^ + E bjX^ 

j=0 



p{xy 



The monic polynomial p{X) is determined by its roots Ai, . . . , Ajv- The polynomial 
q{X) of degree A^ — 1 can be determined from its values at the roots of denominator. 
Therefore, 

Ai, . . . , Xn, qiXi), . . . , g(Aiv); 
are global complex coordinates on this space of maps. 

Let 6 denote a variation of the parameters Ai, . . . , Ajv, ^(Ai), . . . , q{XN), while d 
is a differential of the parameter A. The Atiyah-Hitchin symplectic structure cu is 
defined by the formula 

u = > . AJAfc. 

The corresponding Poisson bracket is specified by canonical relations: 

{?(A„),A,} = 5,"?(A„). (2.1) 

All other brackets vanish 



{A„, Afe} = {q{Xn), q{Xk)} = 0. 
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(2.2) 



The bracket turns the space of maps X{X) : CP^ — > CP^ of degree N with the 
above mentioned properties into a Poisson manifold. 

Consider some function Xq, a point of the Poisson manifold. Fix some points 
A and fi on the sphere away from the poles of Aq. Then the points A and /i 
considered as an argument of X, where X is from a small vicinity of Xq, are 
well defined functions in this vicinity"^. Their values at X are functions of the 
coordinates Ai, . . . , Ajv, ?(Ai), . . . , q{XN)- As it was demonstrated by Faybusovich 
and Gekhtman, [3], the bracket for X{X) and X{iJ,) is given by the formula 

A — H 

Here some miracle occurs. The Poisson bracket for two functions (A and //) on the 
Poisson manifold is given in terms of these functions and their argument {X). In 
[11] for rational functions we gave a direct proof that 2.3 implies 2.1-2.2. We list 
properties of the AH bracket. 

Evidently, 2.3 is skew-symmetric with respect to A and It is also hnear in its 
arguments 

{aX{X) + bX{i2), X{u)} = a{X{\), X{u)} + X{u)}, (2.4) 

where a and b are constants. The Leibnitz rule holds 

{X{X)X{f,),X{u)} = X{X){X{f,), X{u)} + X{f,){X{X), X{u)}. (2.5) 

The bracket satisfies the Jacobi identity 

{X{X), {Xifx), X{u)}} + {X{fx), {X{u), X{X)}} + {Xiu), {X{X), X{fx)}} ^ 0. 

The bracket is compatible with the Cauchy theorem in the following sense. Any 
function X{X) analytic inside some contour C is determined by its the values at 
the contour 

x(0 



c 



C-A 



for A inside C. Whence due to 2.4 the values of the bracket at different points are 

c 

for any fi G CP^. It can be verified directly that the formula on the right in 2.3 
satisfies this compatibility condition. 

The following is particularly useful to us. 



^Instead of X{X) and X{fi) one should write X{X) and /i(Af). 
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Lemma 2.1. The bracket 2.3 is invariant under linear fractional transforma- 
tions 

cX + d' ^ ^ 

where a,b,c,d are constants. 

Proof. Consider a transformation of the form 

X X' ^ J (2.7) 

cX + d ^ ' 

Then 

1 1 {{cX{\) + d)-{cX{^Ji) + d)f 

{cX{X) + df {cX{^Ji) + df A-At 

^ {x'{X)-x'{^,)f ^ 
\ — ji 

To finish the proof we note that two consecutive transformations of the form 2.7 
produce the whole group 2.6. □ 
The Atiyah-Hitchin bracket in coordinate free form 2.3 will appear for a much 
wider class than rational functions. Formula 2.3 itself can be a starting point for 
construction of the bracket on an infinite algebra of complex observables A' (A), A G 
CP^, meromorphicaly dependent on the parameter. Indeed, we start defining the 
bracket for two observables at different points by 2.3 and the extend it to all 
polynomials using 2.4 and 2.5. It can be verified in a long but simple calculation 
that 2.4 and 2.5 imply the Jacobi identity for the AH bracket. It is of interest 
to find all formulas similar to 2.3 with the properties of a Poisson bracket and 
compatible with the Cauchy theorem. 

3. The spectral problem. 

The NLS equation, 
where ip{x,t) is a smooth complex function, is a Hamiltonian system 

with Hamiltonian H — ^ J + dx —energy and the bracket 

r . . f 5A 5B 5A 5B ^ , , 

{A, B] = 2i / ^ — — — - -7^^ — dx. 3.1 

The NLS equation is a compatibility condition for the commutator 
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with^ 



and 



V2 = -y(T3 + Yo 



' 2 




a 
2 



+ 



V 
V' 



We often omit the lower index and write V — V2. 

3.1. The direct spectral transform. The commutator relation produces an 
auxiliary linear problem 

f2{x,X) 

This can be written as an eigenvalue problem for the Dirac operator (see the 
introduction). Let / denote the transposition of the vector / and let /* denote 
the adjoint of the vector /. Let X^(a, b) be a space of vector functions with the 
property 



f'{x,X) = V{x,X)f{x,X), 



f{x,X) 



(3.2) 



/ 



f*{x, X)f{x, X) dx < 00. 



The Weyl solution, [18], 



e{x,y,X) = 



ei{x,y,X) 
e2{x,y,X) 



is the solution of 3.2 which belongs to L'^[y, +00) or L'^{—oo,y], where y is an 
arbitrary point on the line. Due to the theorem of Levitan and Martunov, [7] 
section 8.6, for a continuous potential iIj{x) the Dirac operator is always in the 
limit point case. It means that for A with S^A 7^ there exists one solution from 
L'^[y, +00) and another solution from i^(— 00,7/]. Evidently the Weyl solutions 
are determined up to a multiplicative constant. 

Pick some a G [0,7r). Consider the fundamental system of solutions p^{x,y,X) 
and Ua{x, y, A) of 3.2 normalized by 



Pa{x,y,X)\x=y = 



le 
-ie~ 



"^aix, y, A) \x=y 



(3.3) 



The Weyl solution e{x,y,X) from L'^[y, +00) or L'^{—oo,y] is proportional to a 
linear combination of p^, and Ua- 

e ~ ap„ + bUa, 



^Here and below a denotes the Pauli matrices 



CTl 



ii 1 

1 



(72 = 



-i 

1 



0-3 



1 

-1 
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but the Weyl function X^{y, A) = a/6 is defined uniquely. Evidently, the solution 
X^Pa + Ua belongs to L'^[y, +oo)/L^{—oo, y] for A with S>A 7^ 0. 

Example 3.1. The trivial potential ip{x) = 0. 



_lQ-ia+%{x-y) 



Ur 



{x,y,X) 



The Weyl function X+{y, A) = i if $>A > and -i for S>A < 0. The Weyl function 
X~{y, A) = — i if 5JA > and i for $>A < 0. The only potential which has such 
Weyl functions vanishes identically. 

Lemma 3.2. Any two functions X^{y, A) and X^{y, A) are related by the equa- 
tion 

, A'i^ cosfo; — 3) — sm(a — 3) 
= .r^-^ ^. (3.4) 



sin (a 



cos a 



Proof. First we will obtain the expression for the function X^{y, A) in terms of 
the square integrable solution e{x, y, A). This solution is proportional to X^p^+u^ 
with some constant c. Thus X^p^ + Ua — ce and normalization conditions 3.3 
imply the system 



-X^{y,X)ie-'- + e- 



cei(y,y, A), 
ce2(y,y, A). 



Solving for X^, 



(3.5) 



This identity can be written as 



Xz 



(cie — 626*'^) cos(a — P) — {ciie + 62^6*^) sin(a — f3) 
(eiie~*^ + e2ie^^) cos{a — (3) + (eie~*^ — 626*^) sin(Q; — /3) ' 



Dividing the numerator and denominator by eiie"*^ + 62^6*^ we obtain 3.4. □ 
Therefore, if X^ is known for some value of the parameter a, then it is known for 
all other values of a. 

If f{x, A) is a solution of the auxiliary problem 3.2 corresponding to A, then 
/ = (7i/ is a solution corresponding to A. Formula 3.5 implies, for A with $JA 7^ 0, 



ei(A)e-*" - e2(A)e^' 



e2(A)e---ei(A)e- 



ei(A)ie-^" + e2(A)ie^" e2(A)ie-^" + ei(A)ie^° 

The function X^{y,X) for fixed y maps the upper/lower half-plane into the up- 
per/lower half-plane. 



For a nonzero potential the Weyl function X^{y,X) takes valus in the upper 
hulf-plane when A is in the upper half-plane. It is represented by the integral, [5], 

+00 

X^{y,X)^bX + a+ J 
—00 

with 

+00 

6 > 0, aeR\ 

The spectral measure daa corresponds to a self-adjoint extension of the operator 
D on L'^[y, +00) specified by the boundary condition 

/i(y)e-*" = /2(|/)e^" 

The function X~{y, A) has a similar representation. 

For each point of the complex plane with nonzero imaginary part there are two 
Weyl solutions. To make the Weyl solution a single valued function of a point 
we introduce F = r_|_ U r_, a two sheeted covering of the complex plane (Figure 
1). Each sheet r_|_ or F_ is a copy of the complex plane without the real line. 
Each point of the cover F is a pair Q = (A, ±) where A is a point of the complex 
plane and the sign ± specifies the sheet. We denote by P+ and P_ the infinity 
corresponding to the sheet F+ or F_ (Figure 2). 




Figure 2. Two sheets of the spectral cover. 

Let us introduce two components of the spectral cover 

Tfl = {g e F+, '^X{Q) > 0} U {g e F_, QX{Q) < 0} 

and 

Fl = {g e F+, QX{Q) < 0} U {g e F_, ^X{Q) > 0}. 
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t-X + 1 



daa{t), 



f daa{t) 




Q e r, which is the vector 



Evidently F = Tji U F^. On F we define an involution by the rule 

ea-. (A,±)^(A,T). 

The involution permutes infinities Sa '■ P-\- — > Evidently Fr or F^ are invariant 
under the action of €„. 

We define the Weyl solution e{x,y,Q), a;, |/ G M 
function 

ei{x,y,Q) 
_e2ix, y,Q) 

as a solution of 3.2 with A = \{Q) and the following properties 

• e{x, y, Q) belongs to Li^ly, +00) if Q e F^; 

• e{x,y,Q) belongs to oo,y] if Q £ F^; 
normalized by the condition 



e{x,y,Q) = 



ei{x,y,Q) 



x=y 



= 1. 



(3.6) 



The standard transition matrix M{x,y,X) = (M'^^^M^^)) is 2 x 2 matrix so- 
lution of the auxiliary linear problem 3.2 which satisfies the boundary condition 
M{x, y, X)\x=y — I- The symmetry of the matrix V{x, A) 

aiV(x, A)(7i = V(x, A) 

produces the same relation for the transition matrix 

aiM(x,y,X)ai^ M(x,y,X). (3.7) 

This implies for the columns 

M^'\x,y,X) = a,M'^'\x,y,X), M^''\x,y,X)=a,M^'\x,y,X). (3.8) 

Evidently the Weyl solution e{x,y,Q) is a linear combination of the columns of 
the transition matrix. The normalization condition 3.6 implies 

e{x, y, Q) = M(i) (x, y, A) + X{y, Q)M'^'\x, y, A), (3.9) 

where the function X{y, Q) is called the Weyl function. 

Formulas 3.8, 3.9 and the uniqueness of the Weyl function imply 

1 



X{y:(^aQ) 



X{y,Q) 

Example 3.3. The trivial potential 'ijj{x) = 0. 

For the columns of the transition matrix we have 



(3.10) 



MW(a;,7/,A) =e" 



1 




M(2)(a;,y,A) = e 




1 



We put X{y, Q) = 00 when Q G F+ and X(y, Q) = when Q e F_. 
For t/j{x) 7^ the function X{x, Q) has a pole at P+ and a zero at P_ 
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Arguing as in the proof of Lemma 3.2, for Q G we have 

^(y, Q) - vi^^i ^"'^ A - m- (3.11) 



Note that the transformation 



Z — Z 

estabhshes a 1:1 correspondence between the upper/lower half-plane and the ex- 
terior/interior of the unit circle. This implies, in particular, that the function 
X{y, Q) for Q & with Q=Q > or < takes values in the exterior or interior 
of the unit circle. The shift a — > (3 described on the 2;-plane by formula 3.4 
corresponds to rotation via the angle 2{a — (5) on the 2;'-plane. 
For (5 e we also have 

^(?/,g) = '^5l'^'u ^^'^'"' A = A(Q). (3.12) 

I - [y, A) 

The function X{\,Q) for Q E Tl with '^Q < or '^Q > takes values in the 
exterior or interior of the unit circle. 

We associated to a potential ijj{x) defined on the entire line a pair (F, X). Thus 
we defined a map 

M — > (r,X), (3.13) 
which we call the direct spectral transform. This map is injective. If the function 
X{y,Q) is known for some fixed value of y, formulas 3.11-3.12 allow one to re- 
construct X~{y,X) and X^{y,X) (or equivalently the spectral measure). By the 
theorem of Marchenko, [8], a potential is determined on the corresponding half hne 
by its spectral measure. 

Due to its injective character the direct spectral transform can be inverted 

{T,X) — > M. (3.14) 

This map we call the inverse spectral transform. An effective procedure for the 
inverse spectral transform was constructed by Gelfand and Lcvitan, [4]. 

3.2. The formal series for the Weyl functions. For the function X{y., Q), Q E 
F+, '^Q > 0, and z < y we have 

X(y Q) = ^2(l/,|/,<3) ^ e2{y,z,Q) 
ei{y,y,Q) ei{y,z,Q)' 
This together with 3.2 implies the Ricatti type equation, A = X{Q): 

^ ^xX{y, Q)+^P- ^X^y, Q). (3.15) 

In fact, the function A^d/, Q) satisfies the same equation in the variable y for any 
Q. 
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For a generic potential the Weyl function X{y,Q) has a pole at infinity for 
Q e r_,_ (compare Example 3.3). The formal expansion has the form 

X(y, Q) = a_i(y)A + a^{y) + ^ + . . . , Q e (P+); 

with 

= }j, «o(y) = ^2^, etc. 

The function X{y,Q) has a zero at infinity for Q e r_. It has the power series 
expansion 

XM) = '-^ + '^ + '^ + ..., Qe(P-y, 

where 

The coefficients of the two expansions are obtained by substituting the series into 
the difi^crcntial equation 3.15 and then by matching terms with the same power of 
A. The expansions are connected in accordance with formula 3.10. 

V. A. Marchenko, [9], studied the question of when the formal series have an 
asymptotic character. If the potential has some degree of differentiabihty then the 
formal series have an asymptotic character. According to [9] a converse statement is 
also true. If the formal series converge for large values of A, then the corresponding 
potential is an infinitely differentiable function. In fact, for our purposes it is 
enough to have only the first term of the asymptotic expansion. 

The Weyl functions X^{y, A) and X~[y, A) for fixed A and a are the functions of 
y and the potential. Formula 3.5 and the spectral problem 3.2 imply the Ricatti- 
type equation for X^{y, A): 

^^^^^ = ~{x^Hy+A)+i)-x^{y,x)(^a+^a)+'^{x^ 

where i/ja ~ ijje^'^". In general, the function X[y,Q) is simpler to work with than 

x^iy,^)- 



4. The Poisson bracket. 

4.1. The Atiyah— Hitchin bracket for the Weyl function. The main result 
of this section is Theorem 4.1. This theorem describes the image of the Poisson 
bracket 3.1 under the direct spectral transform 3.13. 

Theorem 4.1. Let y be some fixed point of the line and X{Q) — X{y, Q) and 
X{P)^X{y,P). 
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i. If Q and P belong to Tr, then 

ii. If Q and P belong to Tl, then 

WQ), XiP)} - -2 X -^^(gpAfp)— 
in. If Q belongs to Tr, P belongs to Tl, then 

{x{Q),a:{p)} = o. 

Remark 4.2. In the formulation of Theorem we have to consider various parts of 
the spectral cover separately because it is not simply-connected. 

Formulas 3.11 and 3.12 and the invariance of the AH bracket under linear frac- 
tional transformations (Lemma 2.1) produce the bracket for X^{y, A). We precede 
the proof with two auxiliary results. 

Lemma 4.3. Let the vectors /^(x, A), A) satisfy 

f\x, A) = V{x, \)f{x, A), /♦'(x, A) = V{x, A)/*(x, A), 

and the vectors g^{x,fi), g^{x,fi) satisfy 

g^\x, = V{x, /i)g^(a;, /x), g*'{x, n) = V{x, i^)g*{x, jj). 

The following identity holds: 

d_ 

— A) dx 

Proof. The identity can be verified by differentiation. □ 
Lemma 4.4. If Q E Tr, then the following identities hold: 

^^^"'^^ = -[ei(y,z,g)f, y>z- (4.1) 



SX{z,Q) 



6i){y) 

The derivatives vanish for y < z. 

If Q e Tl, then the following identities hold. 



[e2{y,z,Q)]\ y>z. (4.2) 



Si/jiy) 

The derivatives vanish for y > z. 



^-^^ = hfa,.,0)f, y<z; (4.3) 

SX(z,Q) ^ .le,(y,z,Q)f, y<z. (4.4) 
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The lemma was proved first in [16] . Here we present a simpHfied proof of formulas 
4.1 and 4.2. The proof of other formulas 4.3 and 4.4 is the same. We split the 
proof into small steps. 

Proof. Without loss of generality wc assume that z = 0. 

Step 1. This step prepares gradients. Let M' = 5M be a variation of M{x, 0) in 
response to the variation of ijj{y) and ipiy), < y < x. Then M* = VM* + V'M. 
The solution of this nonhomogenious equation is 



X 



Therefore, 



SM{x, 0) 
5ij{y) 



^M{x,y)(^^ o)M(y,0), 



Step 2. The purpose of this step is to prove the formula 



Sijjiy) 
where 



'mu{y,0) mi2{y,0) 
A B 



(4.5) 



A = m2i{y,0) - X{y,Q)mii{y,0), 
B = -m22{y,0) + X{y,Q)mi2{y,0). 

Consider the auxiliary spectral problem 3.2 on the finite interval [a,b]. The 
solution 

f{x, A) = M^^\x, a, A) + X{b, a, Q)M^'^\x, a, A), 
with some A" (6, a, Q), satisfies the boundary condition A) = /2(&, A) if 



The limit 



A'(6,a,g) = ^^^il^^(6,a,A). 
m22 - mi2 



X{Q,Q) = lim A'(6,0,g) 



exists because the spectral problem is in the hmit-point case. Therefore,^ 

VA'(o,g) = lim vA'(6,o,g). 

6— >+oo 



= S/5tp{y) or 6/5tp{y). 
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To compute the derivative 



we use the formulas of Step 1: 
SX{b, 0,Q) 



mil — ^21 Vmii — Vm2i 
m22 - mil - "^21 

mil - ^121 Vm22 - Vmi2 

m22 — mi2 17122 — mi2 



(b, 0, A) 
(6,0, A). 



(4.6) 



x{b, o,g) 



X 



[mi2 - m22]{b,y)mii{y, 0) 

[mil -m2i]{b,0) 
[m22 - mu] {b, y)mi2(y, 0) 
[m22 -mi2](6, 0) 
X{b, 0, Q) [mi2 - m22] {b, y) 

mii{y,0) mi2(y,0) 



- X{b,0,Q) 



+ 



[mil -m2i](6, 0) ["^22 - "^i2](&,0)_ 
Using the identity M(6, 0) = M(6, y)M{y, 0) and simple algebra, one finds 



_m2i(y, 0) - ^"(6, y, Q)mu{y, 0) -m22{y, 0) + ^"(6, y, Q)mi2{y, 0) 

Now pass to the limit when — > oo. 

5. The purpose of this step is to prove the formula 



5X{0, Q) 
5ijj(y) 



-X{0,Q)X{y,Q) 



^2i(y,0) rn22iy, 0) 
A B 



(4.7) 



To compute the derivative VA'(6, 0, A) given by 4.6 we use the formulas of Step 1 



SXjb, 0, Q) 
5'^{y) 



x{b, o,g) 



- X{b,0,Q) 



[mil - m2i\{b, y)m2i{y, 0) 
[mil -m2i](6,0) 

["'21 - "in]{b. y)ni22{u, 0) 



X 



[m22 -mi2](&, 0) 
X{b, 0,Q)[m,ii-m2i](b,y) 

in2i{y,0) ^22(^,0) 



+ 



[mil - m2i](6, 0) [m22 - mi2](^0)_ 

Using the identity M(6, 0) = M(6, y)M{y, 0), and simple algebra, one finds 

+ A'(6,0,A)A'(6,y,A) 

"^2i(y,0) "7-22 (y,0) 



X 



+ 



_-m2i(y, 0) + ^"(6, y, A)mii(y, 0) m22(y, 0) - X{b, y, A)mi2(y, 0) 
Now pass to the limit when 6 — > 00. 
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Step 4- Consider e{x,y,Q) proportional to e{x,0,Q). Then 



X{y,Q) 



e2{y,y,Q) ^ 62(1/, 0,(3) 

ei{y,y,Q) ei{y,0,Q)' 



(4.8) 



Therefore, 



X{y, Q)[mii(y, 0) + ^"(0, Q)mu{y: 0)] = m2i(y, 0) + ^'(O, Q)m22(?/, 0). 



After simple algebra. 



A ~ B' 



Therefore, using 4.5, 4.8, and 4.9, we obtain 



(4.9) 



SX{0,Q) 
5ip{y) 



= X{0,Q) 

= A'(o,g) 



A 

ei(^,0, A) 

[ei(y,0,Q)f 

m2i(?/, 0)ei(y, 0, Q) - mii(y, 0)e2{y, 0, Q) 



The denominator does not depend on y and can be computed for y — 0, where it 
is equal to — ^"(0, Q). Formula 4.1 is proved. 
Step 5. Using 4.7, 4.8 and 4.9 we obtain 



SXjO, Q) 

5'^{y) 



-X{Q,Q)X{y,Q) 
-X{Q,Q)X{y,Q) 
-A'(0,g) 



m2i(y, 0) + ^'(O, Q)m22(y, 0) 
A 

62(1/, o,Q) 

m2i{y, 0) - X{y, Q)mu{y: 0) 
[e2{y,0,Q)f 



m2i{y, 0)ei(y, 0, Q) - mu{y, 0)62(1/, 0, Q) 



The denominator does not depend on y and can be computed for y = 0, where it 
is equal to —X{0,Q). Formula 4.2 is proved. □ 
Now we are ready to prove the theorem. 
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Proof. We will prove the first formula. All others can be treated the same way. 
By Lemma 4.4, 



{X{Q),X{P)} 



+00 



^ . r 5X^, Q) 8X{y, P) 8X{y, Q) 5X^, P) 
y 

+00 

= 2i y" -eliC, y, g)e?(e, y, P) + e?(e, y, Q)el{^, y, P) d^. 



y 

Using the identity of Lemma 4.3, 

2i 



[e{C,y,QfJe{i,y,P)f\;' 



i{X{P) - X{Q)) 
^ {X{Q) - X{P)f 
\{Q)-\{P) ■ 

The formula is proved. □ 

4.2. Computation of the Poisson bracket for field variables. The inverse 
spectral transform 3.14 maps the AH bracket on Weyl functions given by Theorem 
4.1 to the phase space. The main result of this section is the following 

Theorem 4.5. The AH bracket for the field variables ip{x) and ip{x) is given 
by the formulas^: 

{iiz),ij{y)} = 0, (4.10) 

q{z)My)} = 0, (4.11) 

U^{z),^{y)} = 2iS{z-y). (4.12) 

These identities are an equivalent form of the Poisson bracket 3.1. 

The inverse spectral transform 3.14 is very implicit. However, if we know X{y, Q) 
for all values of the variable y, then the series presented in 3.2 imply 



and 



hin X{Q)X{y,Q)=ti^{y). (4.14) 

V ^-t — 



^The identities are understood in the sence of generalized functions: u{x) = v{x) if for any 
f{x) G Cq° we have J u{x)f{x)dx = J v{x)f{x)dx. 
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The limits are complex conjugate of each other due to 3.10 and 



lim X{Q)X(x,Q) = lim X(eaeaQ)X(x,eaeaQ) ^ lim "^^^"-^^ 



To prove theorem 4.5 using these formulas one needs to compute the bracket 
{X{x, Q), X{y, P)} when x y. Fortunately one needs this formula only asymp- 
totically when one of the points Q or P tends to infinity. 

Lemma 4.6. Suppose Q ^ P± from the imaginary direction. Then 

{X{y, Q), X{x, P)} ~ e-^(«)(--^){A'(x, Q), X{x, P)}. (4.15) 

Proof. The identity 

e(y, z, Q) = M{y, x, X)e{x, z,Q), A = A(Q), 

implies 

^ ^22(y, X, X)X{x, Q) + 77121 (y, X, X) 
mi2{y,x,X)X{x,Q) +mn{y,x,X)' 

If A = ir, r — > ±oo, then 

M(7/, X, A) ~ e-*^/'(^-")"S 

z.e., the transition matrix near infinity behaves like the solution of the free {Yq — 0) 
equation. To see this write the integral equation for the transition matrix 

M(y, X, A) = e-^^/2(^-^)^3 + T (iee"^^/'^^"^^^'lo(0^(C, x, A). 

J X 

In symbohc form M — R + AM, where R is the solution of the free equation, and 
A is the integral operator. Now expand the solution into the Neumann series 

M^R + AR + A^R + .... 

and take the first term of the expansion. This produces the stated asymptotics. □ 

Remark 4.7. We cut the tail of the Neumann series. An estimate for the tail 
for infinitely smooth potentials is given in [17] and for general square integrable 
potentials in [10]. 

The next lemma establishes that the Poisson tensor is real. 

Lemma 4.8. The Poisson brackets for the field variables ■ip{x) and ip{x) are real 



{mM^)} = {i^{y)Mz)}. 
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Proof. We will prove the first identity. The second identity can be proved along 
the same hnes. 

For Q,P & Ffl, using 4.14, we have 



{ij{y),7jj{z)} = -{iij;{y),iij;{z)} 



- Jm^_ {XiQ)Xiy, Q), A(P)A'(^, P)} 



-Jim^ X{Q)X{P){X{y,Q),X{z,P)}. 



We assume that the limit is taken along the imaginary direction. Using Lemma 
4.6, 



hm^ X{Q)\{P)e'~^^^^^'-y'^{X{z, Q), X{z, P)} 



By Theorem 4.1, 



2 



.-. =- hm X{Q)X{P)e^mi^-y^ ^ 

Q,p^p. ^ ' X(Q)-X(P) 

Now we use Q — eaeaQ, P — ea^aP a-nd invariance of under the action of £„. 
Prom formula 3.10, after simple algebra 

. ^ _ ^.^ A(6„g)A(6,P)e-^(-"Q)(--^) ^ {X{z,eaQ) - X{z,eaP)f 
Q,p^p. X\z,eaQ)X\z,eaP) X{eaQ) - X{eaP) 

A(g)A(P)e-^(Q)(-^) {X{z,Q)-X{z,P)r 
Q,™p+ X^{z,Q)X^{z,P) X{Q)-X{P) 

= -Jm^A(g)A(P)e-(^)(-)(^,^} 

- Hm ^ A(g)A(P){^, ^} ^ -^Wl- 

The last line follows from 4.13. □ 

Now we are ready to prove the main result. 
Proof of Theorem 4-5. Due to Lemma 4.8, identities 4.10 and 4.11 are equivalent. 
We compute the bracket 4.10. Let Q, P e Tji and Q ^ P^, P ^ P_ along the 
imaginary direction. Using formula 4.14 and Lemma 4.6 for f{x) e Cq° and y < z 
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we have: 



z z 

J dyf{y){'>P{z),Hy)}^ - J dyf{y){iHz),i'fp{y)} 

-oo — oo 

z 

-^nm^ I dyf{y){X{Q)X{z,Q),X{P)X{y,P)} 

—oo 

z 

- lim X{Q)X{P) f dyf{y){X{z,Q),X{y,P)} 

— OO 

z 

- ^\im^_X{Q)X{P){X{z,Q),X{z,P)} J dyf{y)e 



iX{P)iz-y) 



Let Q,P be such that A(Q) = —ir, A(P) = —2iT. Since X has a zero at P_ we 
have, using Theorem 4.1 when r — > +00, 

Jm^_ X{Q)X{P){X{z, Q), X{z, P)} ^ mKP)^-^^^^^y^^^^P^ = 0(r-^). 
For the integral we have 

z 

—00 

Therefore, 

{^(^),^(l/)} = 0, y<z. 
Using skew symmetry of the bracket and interchanging y and z we have 

{^|;{z),^|;{y)} = 0, y > z. 

Taking the sum of these two formulas, we obtain 4.10. 

Now we compute the bracket 4.12. Let Q,P E Tji and Q — > P+, P ^ P^ 
along the imaginary direction. Then using formulas 4.13, 4.14 and Lemma 4.6 for 
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f{x) e and y < z we have: 

dyf{y){i^{z),^l){y)} = I dyf{y){-i'4j{z),i'4j{y)} 



z 

= lim J dym{^^yKPmy,P)} 

— oo 

z 

—oo 

z 

= lim-^^l^{X{z,Q),X{z,P)} J dyf{y)e-^'^^^^^ 

— oo 

Using Theorem 4.1, 

{X{z, Q), X{z, P)} - 2 X X(gFA(P^ ■ 

Since A' has a pole at P+ and a zero at P_, we have asymptotically 
X{Q)\{P) 



X%z,Q) 



{X{z,Q),X{z,P)}r. 

A(Q)A(P)^2x '^'(^'^^ 



A'2(^,Q) A(g)-A(P) 
Let P = and X{Q) — ir, r — > +oo; then 



Using steepest descent we have 



z 

r.1 



hm / dyf(y)e-^(^-y^^if(z). 



— oo 



Therefore, 

Since the bracket is real, by Lemma 4.8 we have 

{i/j{z),'i/j{y)} ^ -iS{z -y), y < z. 

By the skew symmetry of the bracket and interchanging z and y, 
{jjj(z),'ijj(y)}^i5(z-y), z < y. 
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Taking the sum of 4.16 and 4.17, we obtain 4.12. 



□ 
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